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Abstract 

Surveillance applications require a collection of heterogeneous vehicles to visit a set of targets. In this article, we consider a fundamental 
routing problem that arises in these applications involving two vehicles. Specifically, we consider a routing problem where there are two 
heterogeneous vehicles that start from distinct initial locations, and a set of targets. The objective is to find a tour for each vehicle such that 
each of the targets is visited at least once by a vehicle and the sum of the distances traveled by the vehicles is a minimum. We present a 
primal-dual algorithm for a variant of this routing problem that provides an approximation ratio of 2. 
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INTRODUCTION 

Heterogeneous unmanned vehicles are commonly used in surveillance applications for monitoring and tracking a set of 
targets. For example, in the Cooperative Operations in Urban Terrain project pQ at the Air Force Research Laboratory, 
a team of unmanned vehicles are required to monitor a set of targets and send information/video about the targets to 
the ground station controlled by an human operator. The human operator may further add new locations of potential 
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£ .targets or task the vehicles to revisit the targets at different angles. Once the human operator enters his/her input 
' through the human-machine interface, the central computer associated with the interface has few minutes to determine 
,— ' \ the motion plans for each of the vehicles. A fundamental subproblem that has to be solved by this computer is the 
, problem of finding a tour for each vehicle so that each target is visited at least once by some vehicle and an objective that 
• • ■ depends on the distances traveled by the vehicles is a minimum. A common objective that is used for these applications 
1 is the sum of the total distances traveled by all the vehicles. If there is only one vehicle, this routing problem is referred 
5— i ' to as the Traveling Salesman Problem (TSP) in the literature. If there are multiple vehicles that (possibly) start from 
different initial locations or depots, then this routing problem is referred to as the Multiple Depot, TSP. Once the 
routing problem is solved and the tours have been determined, a nominal trajectory can be specified for each vehicle to 
include other kinematic constraints of the vehicles using the results in [2] , [3] . 

A multiple depot, TSP is a generalization of the single TSP and is NP-Hard. This routing problem is further 
complicated if the vehicles involved are heterogeneous. In this article, vehicles are considered to be heterogeneous if 
the distance to travel between any two targets depend on the type of the vehicle used. In the context of unmanned 
applications, as a multiple depot heterogeneous TSP is generally a subproblem that needs to be solved, we are interested 
in developing fast algorithms that produce approximate solutions than find optimal solutions that may be relatively 
difficult to solve. In addition, the cost of traveling between any two targets may not be accurately known because 
the exact locations of the targets may not be accurately known, or because the wind conditions used to determine the 
travel cost could change. For these reasons, the main focus of this article is to develop approximation algorithms for 
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heterogeneous TSPs. An approximation algorithm for a problem is an algorithm that runs in polynomial time and 
produces a solution whose cost is at most a given factor away from the optimal cost for every instance of the problem. 

The objective of this article is to develop a primal-dual algorithm for a two depot, heterogeneous TSP (2DHTSP). In 
addition to assuming that the costs satisfy the triangle inequality for each vehicle, we consider a variant of the problem 
where the cost of traveling between any two targets for the first vehicle is at most equal to the cost of traveling between 
the same targets for the second vehicle. Using these assumptions, we show that the developed primal-dual algorithm 
has an approximation ratio of 2. We are motivated to address this variant of the 2DHTSP due to the following reasons: 

• The 2DHTSP considered in this article is one the simplest cases of the general multiple depot, heterogeneous TSP. 
The objective of this work is to first develop good algorithms that can handle these simpler cases efficiently. 

• Consider a scenario where each of the vehicles is modeled as a ground robot that can move both forwards and 
backwards with a constraint on its minimum turning radius. If the minimum turning radius of the first vehicle is at 
most equal to the minimum turning radius of the second vehicle, it follows that the optimal distance required to travel 
between any two targets for the first vehicle will be at most equal to the optimal distance required for the second vehicle. 
Therefore, the problem addressed in this article is a useful variant to address. 

• The 2DHTSP considered in this article is a generalization of a 2 depot, homogeneous TSP where there are additional 
vehicle-target constraints which require one of the vehicles to necessarily visit a given subset of targets in addition 
to visiting any common target available for both the vehicles. This variant of 2 depot, homogeneous TSP arises in 
applications where the distance to travel between the targets are identical for both the vehicles, but one of the vehicles 
carry sensors that require the vehicle to visit a subset of targets compulsorily. 

• For some missions involving identical vehicles, it is sometimes necessary to minimize the maximum cost incurred by 
any of the vehicles. This problem is referred to as the min-max, multiple depot, homogeneous TSP in the literature. 
If there are only two vehicles involved, an application of Lagrangian duality transforms this min-max problem to the 
variant of heterogeneous TSP considered in this article. 

Without the assumptions on the costs of the two vehicles, the 2DHTSP is a generalization of the standard variant 
of the prize collecting TSP considered by Goemans and Williamson in \IU\. In this variant, each target essentially has 
a penalty associated with it. The objective of the prize collecting TSP is to find a tour for the vehicle that starts and 
ends at the depot such that the cost of the tour plus the sum of the penalties of each target not present in the tour is a 
minimum. For any two vertices i and j, if 7Tj, ttj denote the penalties of i and j respectively, then one can pose the prize 
collecting TSP as a 2DHTSP by setting the cost of traveling the edge joining vertices i and j for the second vehicle to 
be equal to ZiizZEa _ Essentially, by choosing the penalty variable corresponding to the second depot to be equal to 0, one 
can deduce that the travel cost for the second vehicle is actually equal to the sum of the penalties of the targets not 
present in the tour of the first vehicle. Even though there are no penalties explicitly mentioned in the 2DHTSP, the 
tour cost for the second vehicle which essentially account for targets not visited by the first vehicle act as penalties. For 
these reasons, the primal-dual algorithm presented in this article is related to the primal-dual algorithm available for 
the prize-collecting TSP in [10| . Essentially, our algorithm is based on the well known moat growing procedure proposed 
by Goemans and Williamson in [TO] . 

Most of the work in the literature related to approximation algorithms for multiple depot, TSPs deal with identical 
vehicles. For example, when the costs satisfy the triangle inequality, there are several approximation algorithms for the 
multiple depot, homogeneous TSP in [3]- [6]. Recently, a 3— approximation algorithm was presented for a two depot, 
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heterogeneous TSP in 7 . This algorithm partitions the targets by solving a linear programming relaxation and then 
uses Christofidcs algorithm [8] to find a sequence of targets for each vehicle. 

The 2-approximation algorithms available in the literature for the multiple depot, TSP generally follow a two-step 
procedure. In the first step, a constrained forest problem which is generally a relaxation of the multiple depot, TSP is 
solved optimally. In the second step, an Eulerian graph is found for each vehicle based on the constrained forest. From 
these Eulerian graphs, a tour can be found for each vehicle by short-cutting any target already visited by a vehicle. 
In this article, we follow a similar procedure where we first find a heterogeneous spanning forest using a primal-dual 
algorithm by solving a relaxation of the 2DHTSP. Then, the edges in the heterogeneous spanning forest are doubled to 
obtain an Eulerian graph for each vehicle. Given these Eulerian graphs, one can [3] always find a tour for each vehicle 
that visits each of the targets exactly once. The crux of this procedure depends on finding a good heterogeneous spanning 
forest. Using a primal-dual algorithm, we find a heterogeneous spanning forest whose cost is at most equal to the optimal 
cost of the 2DHTSP in polynomial time. Hence, it follows that the approximation ratio of the proposed procedure is 2. 



Let D = {di, c^} denote the two depots (initial locations) corresponding to the first and the second vehicle respectively. 
Let T be the set of targets to be visited by both the vehicles. As there are two heterogeneous vehicles, we use two edges 
to connect any two targets such that each edge is used by one of the vehicles. Let V\ := T[J{di} be the set of vertices 
corresponding to the first vehicle and Vi :— TU{d%} be the set of vertices corresponding to the second vehicle. Let E\ 
denote the set of all the edges used by the first vehicle that join any two vertices in V\. Similarly, let Ei represent the 
set of all the edges used by the second vehicle that join any two vertices in V2. Note that E\f]E2 =0. Given any two 
vertices we use {i,j}i £ E\ to represent the edge used by the first vehicle to travel between i and j. Similarly, we 
use {i,j}2 £ E2 to represent the edge used by the second vehicle to travel between i and j. 

Let the cost of traversing an edge e £ E\ for the first vehicle be denoted by cost], . Similarly, let the cost of traversing 
an edge e £ E% for the second vehicle be denoted by costj,. We will assume that it is always cheaper to travel between 
any two targets using the first vehicle as compared to using the second vehicle, i.e., for any two targets i,j £ T, 
cost 1 ^ < cost% We also assume that the costs satisfy the triangle inequality for both the vehicles, i.e., for 
m = 1, 2, for every i, j, k £ V m , e\ := {i, j} m , e 2 ■= {j, k} m , e 3 := {i, k} m , cost™ + cost™ > cost™. 

A tour for the first vehicle starts from its depot d\, visits a set of targets in a sequence and finally returns to d\. A 
tour for the second vehicle starts from its depot c?2, visits a set of targets in a sequence and finally returns to c?2- Note 
that a tour for the first vehicle uses edges only from E\ and a tour for the second vehicle uses edges only from E2. The 
objective of the 2DHTSP is to find a tour for each vehicle such that each target is visited exactly once by some vehicle 
and the sum of the cost of the edges traveled by both the vehicles is a minimum. 



I. Problem Statement 



II. Problem Formulation 



The 2DHTSP can be formulated as an integer linear program as follows: 




(1) 
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x e + 2 z v >2 VS C T, (2) 

e£«i(S) TDUDS 

y^ 2 E z <j vsct, (3) 

ee<5 2 (S) TDUDS 

^ i e +2 ^ Z[/ = 2 V«eT, (4) 

e£<5i(«) T3tO{u} 

^ y e = 2 £ 2( / v«er, (5) 

e£S 2 (u) TDUD{u] 

zu < 1, (6) 

UCT 

x e G {0, 1} for e = {u, w}i Vit, u£T, 

y e £ {0, 1} for e = {u, v} 2 , Vit, u G T, 
a; e G {0,1,2} for e = {di,u}x Vu £ T, 
t/ e €{0,l,2}fore = {(i 2) ?;} 2 ,V«eT, 

zc; G {0, 1} VU C T. (7) 



In the above formulation, x e is the integer variable that represents whether edge e G E\ is present in the tour 
corresponding to the first vehicle. For any edge e joining two targets, x e can take values only in the set {0, 1}; x e = 1 
if e is present in the tour of the first vehicle and x e = otherwise. We have formulated the problem such that any of 
the vehicles can choose to visit just one target if required. Therefore, for any edge e joining the depot d\ and a target 
v G T, x e can choose any of the values in the set {0, 1, 2}. Similarly, y e is the integer variable that represents whether 
edge e G i?2 is present in the tour corresponding to the second vehicle, zjj is a binary variable that determines the 
partition of targets connected to the first and the second depot; zjj is equal to 1 if each target in U C T is connected 
to the second depot using edges from E2 and each target in T \ U is connected to the first depot using edges from E\ . 
Equation ([6]) ensures that zjj can be equal to 1 for at most one subset of T. 

In equations (|2I3I) . Si(S) (for i = 1,2) denotes the subset of all the edges of Ei with one end in S and an other end in 
Vi\S. Si(S) is also referred to as the cut set of S corresponding to the i th vehicle. Equation @ states that for any SCT 
at least two edges must be chosen from S\ (S) if there is at least one vertex in S that is not required to be connected to 
the second depot, i.e., J2e£S 1 (S) > 2 if J2tduds z u = Similarly, equation ([3]) states that for any S C T, at least 
two edges must be chosen from 8%(S) if all the vertices in S are required to be visited by the second vehicle. Equations 
(|4I5|) state that the degree of each of the targets must be equal to 2. In the above formulation, we have not included 
the degree constraints on the depots for the following reason: as the costs satisfy the triangle inequality, one can always 
shortcut the edges incident on the depots in any solution to ensure that the degree of a used depot is 2. Now, consider a 
Linear Programming (LP) relaxation of the problem where the constraints in equations (|4I5I6I7[) are relaxed as follows: 

Ci p = min cost\ x e + cost 2 e y e (8) 



^ i e + 2 ^ z y >2 VSCT, 

ee<5i(S) TDC/DS 

J] J/e>2 ^ VSCT, 

eG<5 2 (S) TDUDS 

x e > Ve e Ei, y e > Ve G E 2 , 
z v > VU C T. 

A dual of the above LP relaxation can be formulated as follows: 

C dual = max 2 ^ ^(5) (9) 

SCT 

E Fi(S)<cos^ Vee^i, (10) 

S:e£5i(S) 

E < cos^ Ve G £ 2 , (11) 

S:eeS 2 (S) 

J2 ^ E VC/ C T, (12) 

SCU scu 

Yx(S),Y 2 {S) >0 VSCT. (13) 

We use the above dual problem to find a Heterogeneous Spanning Forest (HSF). A HSF is a collection of two trees 
where the first tree spans a subset of targets and d\ using edges only from E\ , and the second tree connects the remaining 
set of targets to c?2 using edges only from E 2 . In the next section, we present a primal-dual algorithm that finds a HSF. 
We later show that the cost of this HSF is at most equal to the optimal cost of the above dual. This leads to a 
2-approximation algorithm for the 2DHTSP. 

III. A Primal-Dual Algorithm for finding a Heterogeneous Spanning Forest 

The initialization, the main steps and the final pruning step of the primal-dual algorithm are presented in Algorithms 
[1] , [2] and [3] In the following discussion, we first aim to highlight the important features of the algorithm and give an 
outline of the important steps. The proofs and the correctness of the algorithm will be discussed later. At any iteration, 
the algorithm maintains two sets of connected components corresponding to the two vehicles. In the algorithm, C\ 
denotes the set of connected components corresponding to the first vehicle. Similarly, Ci represents the set of connected 
components corresponding to the second vehicle. For i = 1,2, we use Fi to denote the set of edges connecting the 
vertices in Ci. Initially, both C\ and C 2 consist of components where each vertex is in its own connected component, 
i.e., C\ = {{v} : v G Vi} and C 2 — {{v} : v € V 2 }. That is, both F\ and F 2 are empty. All the components 
in C\ and C 2 are initially active except the components that contain the depots. Also, each vertex in V\ is initially 
unmarked. At the beginning of each iteration, for any VC € Ci, the internal variable w(C) keeps track of J2scc ^i^)' 
i.e, w(C) — Y^scc^iiS)- Similarly, VC G Ci, Bound(C) keeps track of X^scc Y 2 (S). Essentially, variables w{C) and 
Bound(C) are used to enforce the constraints in (fT2"j). At the start of the first iteration, both w(C) and Bound(C) are 
set to zero, i.e., w{C) := and Bound(C) :— VC G C\ {Refer to the initialization steps in algorithm^. 
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As the primal-dual algorithm follows a greedy procedure, and since the cost of traveling between any two targets is 
cheaper for the first vehicle (i.e., cost 1 ^ < cost 2 ^ ^ G T), the algorithm prefers to connect any two targets using 
edges from E\ before using any edges from Ei. Therefore, as the algorithm progresses, the components in C\ tend to 
merge and grow bigger than their corresponding components in Ci. In this article, we refer to the components in C\ as 
parents and the components in Ci as their children. For components C\ G C\ and Ci G Ci, we define C\ as the parent 
of C2 and C2 as a child of C\ if Ci C C\ and di C%. For any component C\ G C\, we use Children(Ci) to denote 
all the children of C\ present in Ci. For any component Ci G Ci,di ^ Ci, we use Parent(Ci) to denote the parent 
of Ci present in Ci. According to the definition, if Ci contains the depot di, Ci doesn't have a parent; however, to 
simplify the presentation, we let Parent(Ci) be an empty set if Ci contains di. At the start of the algorithm, for any 
target v G T, Children({v}) is assigned to be equal to {v} and Parent({v}) is assigned to be equal to {v}. Also, the 
components that consist of just the depots neither have a parent or a child (Refer to the initialization steps in algorithm 

In each iteration of the algorithm, the dual variable corresponding to each of the active components in C\ and Ci are 
increased as much as possible by the same amount until one of the constraints stated in (|10M12[) become tight (Refer to 
lines 2-5 of the algorithm^. If a constraint in (|10[) becomes tight for some edge e G E±, Fx is augmented with this new 
edge and the two components (say C\ X ,C\ V in C\) connected by e are merged to form a single connected component. 
The children of each of the two components C\ x ,C\ y now together become the children of the resulting component 
CixlJCiy. The resulting component becomes inactive if it contains the depot d\, otherwise, it is active. In the case 
when the resulting component becomes inactive, all the children of the resulting component also become inactive (Refer 
to lines 17-27 of the algorithm^. 

Algorithm 1 Primal-dual algorithm: Initialization 
Fx 4- 0; F 2 «- 

& <- {{v} : v G Vi}\ Ci <- {{v} :veV 2 } 

for v G V\ do 

Unmark v 

p(v) <- 

w({v}) <- 

Bound({v}) <— 

If v = di, then Children({v}) <s— 0, else Children({v}) <— {v} 

If v = di, then active\({v}) = 0, else active\({v}) = 1 

end for 

for v G V2 do 

q(v) <- 

H v = d 2 , then Parent({v}) ^— 0, else Parent({v}) <G- {v} 
If v = c?2, then actwe2({f}) = 0, else actii^^f}) = 1 
end for 



Similarly, if one of constraints in becomes tight for some edge e G E2, F2 is augmented with this new edge and 
the two components (say Ci x ,Ci y in C2) connected by e are merged to form a single connected component (Refer to 
lines 29-4-2 of the algorithm^. The resulting component becomes inactive if it contains the depot di; otherwise, it is 
active. In the case when the resulting component is active, the parent of either Ci x or Ci y is assigned as the parent of 
the resulting component (It turns out that due to our assumptions on the costs, when the algorithm enters this part 
of the implementation, both Ci x and C'i y must be active and must be the children of the same parent; we will show 
this result later in the article). In the case when the resulting component becomes inactive, and say C 2x was the active 
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Algorithm 2 : Primal-dual algorithm - Main steps 

1: while 3C E C\ such that active\{C) = 1 do 

2: Find edge e\ — {i,j}i E E\ with i G C\ x ,j E Ci^ where Ci x ,Ci y £ C^Ci^ ^ Ci y that minimizes 

_ ( cost l 1 -P(i)-P(j)) 
1 activei (Cix)Eactivei (Ciy) 

3: Find edge e 2 = {i,j}2 € ^2 with i G C 2x ,j S C^ where C 2x ,C 2y S C2,C2a; 7^ C^j, that minimizes 

_ (cos*e 2 -q(')-q(j)) 

^2 active-2 (C2x)+active2 (C2y) 

4: Find a maximal £ C Ci such that for all C G £ the following properties hold: active\{C) = 1 and Children(C) is 
an empty set. If £ is empty, then e 3 E- 00, else find C G £ that minimizes £3 = Bound(C) — w(C) 

5: e min = min(e-i, £2,2-3) 

6: for C E C\ do 

7: w(C) 4- w(C) + e min x activei(C) 

8: For all v G C, 4— p(v) + e m in x active\{C) 

9: Bound(C) E- Bound{C) + J2deChiidren(c) £mm x active 2 (C) 

10: end for 

11: for C G C 2 do 

12: For all v G C, q(w) 4— g(w) + e m i„ x active 2 {C) 

13: end for 

14: switch £ mm 

15: //Comment: If more than one value in {£1,62,63} is equal to e m i n , then give priority first to Case 61, 
then to Case £ 2 and finally to Case £3 

16: Case 61: 

17: Fx^-flUM 

18: Cl^CiUICl^UClJ-Clz-Cia 
19: wCCizUCly) 4-w(Ci x )+w(Cij / ) 

20: Children(Ci x [j C\ y ) E- Children(C lx ) [j Children(Ci y ) 
21: For all C E Children{C\ x [j Ci y ), Parent(C) E- C^ 1J Ci v 
22: Bound{C\ x 1J Ciy) <— Bound(Ci x ) + Bound(Ci y ) 
23: if di G Ci^lJCiy; then 
24: activei(Ci x [J C\ y ) = 

25: active 2 (C) = for all C G Children(C\ x \] C ly ) 

26: else acizi;ei(Ci a ; 1J C\ y ) = 1 

27: end 

28: Case e 2 : 

29: F 2 -E- F 2 {J{e 2 } 

30: C 2 <— C 2 lJ{C2a; U C 2 y} — C 2x — C 2y 

31: if d 2 E C 2x U C 2y then 

32: active 2 (C 2x [j C 2y ) E- 
33: Parent(C 2x [j C 2y ) E- 

34: Let C G {C 22: ,C 2y } such that d 2 ^ C; Children(Parent(C)) E- Children(Parent(C)) - C 
35: else if active 2 (C 2x ) — 1 and active 2 (C 2y ) = 1 
36: active 2 (C 2x \J C 2y ) E- 1 

37: C temp 4- Parent(C 2x ) 

38: //Comment: Parent(C 2x ) must be the same as Parent(C 2y ) 
39: Parent(C 2x [j C 2y ) E- C temp 

40: Children(C temp ) E- Children(C temp ) \J{C 2x \J C 2y } - C 2;E - C 22/ 
41: end if 
42: end if 

43: Case £3: 

44: activei(C) E- 

45: Mark all the unlabeled vertices of C with label C 

46: end switch 

47: end while 
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Algorithm 3 : Primal-dual algorithm - Pruning step 

1: F[ is obtained from Fi by removing as many edges as possible from _F\ so that the following properties hold: 1) All 
the unmarked vertices of V\ are connected to the first depot d\\ 2) If any vertex with label C is connected to the 
depot d\ , then any other vertex with a label C 3 C is also connected to the depot d\ . 

2: F 2 is obtained from F 2 by removing as many edges as possible from F 2 such that any target not spanned by F[ is 
spanned by F 2 . 



component during the iteration which did not contain the depot, the parent of Cix looses Cix as its child. 

If one of the third set of constraints in (fT2"|) become tight for some component C in C\ , then the algorithm deactivates 
C and marks each of the unmarked vertices in the component with C {Refer to lines 44~45 of the algorithm^. The 
algorithm terminates when all the components in C\ become inactive. After termination, the algorithm makes one final 
pass at all the edges (refer to algorithm^ and removes any edge that is not required to be in the HSF. Basically, during 
the final step of the primal dual algorithm, any unnecessary edges in F\ and Fq, are pruned further to find a tree for 
each of the vehicles. Specifically, the tree corresponding the the first vehicle is obtained from F\ by removing as many 
edges as possible from F\ so that the following properties hold: 1) All the unmarked vertices of V\ are connected to the 
first depot d\; 2) If any vertex with label C is connected to the depot d\, then any other vertex with a label C" 3 C is 
also connected to the depot d\. The tree corresponding the second vehicle is obtained from F% by removing as many 
edges as possible from F 2 such that any target not spanned by F[ is connected to d 2 in F 2 . 

Since the sum of the number of components in C\ , the number of active components in C\ and the number of com- 
ponents in C2 decrease at least by one during each iteration, the primal-dual algorithm must terminate after at most 
3|T| + 2 iterations. Using the techniques given in [10], this primal-dual algorithm can be implemented in \T\ 2 log \T\ steps. 

Remark: Case £2 is never chosen by the algorithm if one of the components in Cqxi Ci v is inactive and at the same 
time, the depot d 2 is not present in either Cq x or C^y For this reason, the if loop in lines 31-42 does not consider this 
case at all. The proof of this statement is shown in lemma [TJ 

Remark: In lines 36-40 of the main algorithm, an active component G C2 merges with another active component 
Ciy G C2. We claim that the algorithm will enter these lines only if both C2X and C2 V have the same parent. We show 
this in lemma^ 

Remark: At the termination of the algorithm, F[ will connect a subset of targets using edges from E\ to the depot 
d\. F 2 will connect each of the targets not spanned by F[ to depot using edges from E2. The proof of this statement 
is shown in lemma® 

In the following subsection, we first discuss some properties of this primal-dual algorithm before proving its approxi- 
mation ratio. 



9 



A. Properties of the primal-dual algorithm 

Observation 1: At the start of an iteration of the algorithm, for any two disjoint components C\ x , C% y <E C\, consider 
the constraint in (|10p corresponding to the edge e = {u,v}i that could potentially connect vertex u in C\ x to vertex 
v in C\ y : ^2s-e&Si (S) 0^) — cos ^l- Since e has not yet been added to Ft, this constraint can be re- written as 
Ssnes^ L (^) + Yls-veS^iiS) < cost\. To simplify this constraint further, let us define p(v) = Yls-vesYi(S) Vw S Vi. 
Now, the constraint can be formulated as p(u) + p(v) < cost],. 

We can interpret as the total price vertex u is willing to pay to buy any edge that is incident on vertex v. There- 
fore, to add an edge {it, v}i during the iteration, both the vertices u and v must increase the sum of the prizes they are 
willing to pay by cost] —p{u) —p(v). Hence, each of the dual variables of the active components have to be increased by 
an amount given by act j„ei(Ci )+activei(Ci ) * n or< ^ er to make the constraint, p(u) +p(v) < cost], tight. Hence, in step 2 
of the algorithm [21 we basically find the minimum amount by which each of the dual variables of the active components 
in C\ have to be increased so that none of the constraints are violated and at least one of the constraints in (|10[) just 
become tight. Similarly, in step 3 of the algorithm [21 we find the minimum amount by which each of the dual variables 
of the active components in C2 have to be increased so that none of the constraints are violated and at least one of the 
constraints in (fTTj) just become tight. 

Observation 2: For any active component G\ G Ci, as long as C\ has at least one child, the constraint associated with 
C\ in p2[) never gets tight. In the algorithm, the variable Bound(Ci) keeps track of the amount by which J^scd ^(S) 
is increased during each iteration. Once an active component C± looses all its children, Bound(C\) specifies the maxi- 
mum value that can be attained by io(Ci) = X^scci ^i(^)- 

Observation 3: An active component in C\ can get deactivated during an iteration due to one of the following reasons: 
1) the component merges with another component in C\ that consists of the depot di, 2) the dual variable of the active 
component cannot be increased more than a certain value as it would violate its associated constraint in (j 12|) . On the 
other hand, an active component in Ci can get deactivated during an iteration due to one of the following reasons: 1) 
the component merges with another component in C2 that consists of the depot d<i, 2) the parent of the component gets 
deactivated, i.e., the parent of the component gets connected to depot d\. 

Observation 4-' If an active component C± G C± merges with another component containing d\, then the merged 
component and all the children of the merged component (which includes the children of C\ ) also get deactivated. For 
this reason, each of the children of any parent containing d\ is always inactive during any iteration of the algorithm. 
However, if a component C2 € C2 merges with another component containing c?2 and gets deactivated, then the parent of 
C2 looses C2 as its child but can still be active. Moreover, the dual variable corresponding to parent of C2 can continue 
to increase until one of the constraints in ([T0| or (fT2|) becomes tight. 

Observation 5: In each iteration, for any target v G T, p(v) is incremented by e m in if v belongs to a component in C\ 
that is active; else p(v) does not change. Similarly, for any target v £ T, q(v) is incremented by e m i n if v belongs to a 
component in C2 that is active; else, q(v) does not change. 
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Lemma 1: Consider any target u G T. At the start of the k th iteration, let C\ (u) denote the component in C\ 
containing u, and Cf (it) represent the component in C2 containing u. Then, each of the following statements are true: 

1. Consider a case when C\ (it) is a parent of Cf (it). If C'l(u) is inactive because its parent C\ (it) contains d\, then 
C\ (it) can never merge with any other component during the k th iteration. 

2. For any two distinct targets it and v, if components Cf (11) € C2 and C^v) € C2 are active, then C|(it) can merge 
with (i>) during the k th iteration using an edge {u,v}2, only if both C\ (it) and C\ (v) have the same parent. 

3. Consider a case when C\(u) is active and is a child to an active C\ (11). Then, at the end of the k th iteration (or 
at the start of the k + 1 th iteration), the parent C^ +1 (u) does not have any child which contains target u only if Cf (u) 
merges with another component in C2 that consists of the depot d,2 ■ 

4. If active\{Ci{u)) > activeziC^ (it)) , then at the end of the k th iteration (or, at the start of the next iteration), 
activei(C% +1 (u)) > active 2 {C^ +1 (it)) . 

Proof: Let us prove this lemma by induction. During the first iteration, there is no component in C2 that has 
a parent consisting of d\. Therefore, lemma 1.1 is true by default. At the start of first iteration, each of the active 
components in both C\ and C2 contain only one target. Also, for any u G T, the component {11} G C\ is the parent of 
{11} G C2. Since the cost of the edge joining any two targets for the first vehicle is at most equal to the cost of the edge 
joining the same targets for the second vehicle, the algorithm will always prefer to edge {it, v}i as compared to {it, v}2- 
Therefore, no edge corresponding to the second vehicle joining two targets will be added during the first iteration of the 
algorithm. Hence, lemma 1.2 is true by default for k = 1. The only way the parent {u} G C\ looses {u} G C2 as the 
child is if {u} G C2 merges with the component {0I2} during the first iteration. In this case, a new inactive component 
{it, (^2} G C2 is formed and {u} G C% looses {11} G C2 as the child. Hence lemma 1.3 is correct for k = 1. At the start of the 
first iteration, each of the components which contain one target is active. Therefore, activei{C\(u)) > active2(C2(u)). 
Also, during the first iteration, one of the following cases can happen: (a) an edge is added to F\ joining two targets or 
a target and d±, (b) an edge is added to F2 joining a target and c?2- In either case, one can verify that the components 
are merged such that active\{Ci{u)) > acbiveziC^u)) for all u G T. Hence lemma 1.4 is correct for k = 1. 

Now, let us assume that all the statements in lemma 1.1-1.4 are true for the I th iteration for any / = 1, • • • , k — 1. 
From lemma 1.4, this assumption implies that active i(C[(u)) > active2(C 2 (u)) for any I = 1, • • • , k. We then show that 
all these statements are also true for the k th iteration. 

Proof of lemma 1.1: In this case C\ (it) is a parent of C|(it), and C\(u) is inactive because its parent C\ (it) contains 
d\. Since C\(%i) contains di, from observation^ it follows that all the children of C\ (u) are inactive. Therefore, if C\ (u) 
which is inactive has to merge with some other component C\ (v) corresponding to target v, then Parent(C2(u)) 
Parent^C^iv)) and C^iv) must be active. Since activ ei(C[(u)) > active2(C2(u)) for any I = 1, ■ • • , fc, from observation 
[U it follows that p(u) > q(u) at the start of the k th iteration. Similarly, for target v, we must also have p(v ) > q(v) at 
the start of the k th iteration. Since, cost) < cost 2 , 1 , we now have 

COsf {«,u}i -P( U ) -P( V ) 

1 activei(C^(u)) + activ ei(C^(v)) 

cost lu,v}2 ~ ~ q ( v ^ 
' active2(C2 {u)) + acttwe2(C|(u)) 



,.. S-ZT- ^7^ = ^ (W) 



As a result, the algorithm will prefer to add edge {u, v}i as compared to adding {u, v}2- But, once {u,v}i is added, 
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the component C^iv) will become a child of (u) [j C\ (v) and as a result will be deactivated. Therefore, C^u) will 
remain inactive and will never merge with any other component during the k th iteration. 

Proof of lemma 1.2: For each target u £ T, if activei(C[(u)) > active2{C l 2 {u)) for any I = 1 • • • , fc, then from 
observation [5j it follows that p(u) > q(u) at the start of the k th iteration. Let C 2 (u) £ C2 and Cf (t>) G C2 be active and 
disjoint components, and let Parent(C 2 (it)) ^ ParentiC^iy)). As cosi| u K i < cosi| u ^ for any u,v £ T, again using 
equation I14[ we can conclude that e\ < £2- If £1 < £2, edge {w,i;}2 cannot be added during the iteration. If E\ = £2, 
as mentioned in line 15 of algorithm [2 preference is given to adding edge {u,v}i as opposed to adding edge {it,u}2- 
Therefore, in either case, edge {it, v}2 cannot be added during the iteration. Hence, two active components in C2 can 
merge only if they have the same parent. 

Proof of lemma 1.3: Due to the manner in which the components are deactivated in the algorithm {Refer to line 
34 in algorithm^ and observation^, if Cf (it) is active and is a child to C±(u), at the end of the k th iteration, the 
parent C^ +1 (u) does not have any child which contains target u only if C|(u) merges with another component in C 2 
that consists of the depot 0I2 ■ 

Proof of lemma l.^: 

If C 2 (u) is inactive, from observation [3j either C 2 {u) must contain the depot 0I2 or its parent C\ (it) must contain 
the depot d\. 

• If C|(it) already contains c^, then C*2 +1 (u) must also be inactive. Therefore, activei(Ci +1 (u)) > active2(C2 +1 (u)) 
0. 

• If C\ (u) is inactive because its parent C\ (u) contains di, then we have already shown that C^u) can never merge 
with any other component during the k th iteration. Therefore, activei(Ci +1 (u)) > active2(C2 +1 («))■ 

If C^u) is active, then active\(Ci(u)) > active2{C\ (u)) implies that C\{u) is also active. In addition, since the 
lemma is true for any I = 1, • • • , k — 1, from lemma 1.3 it follows that C\ (u) is a parent of C^iu). Since the component, 
C\ (u), has at least one active child in C|(u), C\ (u) can never become inactive due to its associated constraint in fj 12|) 
during the k th iteration. The only way (u) can lead to an inactive Cf +1 (u) is if (u) merges with another component 
containing di during the iteration in which case all the children of C± (u) including C|(u) also get deactivated. Therefore, 
activ ei (C* +1 (u)) > active 2 (C^ +1 (u)). 

m 

Corollary 1: An edge joining two targets u,v £ T can be added to F2 only if both u and v are already connected in 
JPi, i.e., both u and v belong to the same component in C±. 

Proof: We use the same notation we used in lemma [TJ If target u is not connected to v at the start of the k 
iteration in both (Vi,.Fi) and (V-^-Fz), then using the results from lemma [T] and equation [TU £1 < £2. Therefore, the 
algorithm always prefers to connect targets u,v first using {u, v}i as compared to using {u,v}2- Hence proved. 

■ 

Lemma 2: The algorithm produces a feasible, heterogeneous spanning forest, i.e., the trees specified by the collection 
of edges in F[ and F' 2 connect each of the targets to one of the depots. Any vertex spanned by the edges in F[ is not 
spanned by the edges in F' 2 and vice versa. 

Proof: The algorithm terminates when all the sets of C\ become inactive. This is only possible if each of the targets 
in T is either connected to d\ or c?2. Note that F[ is formed from F\ such that each of the unmarked vertices remain 



12 

connected to d\. The only vertices not spanned by F[ are some of the marked vertices. These vertices were marked 
because the components in C\ that span these vertices were deactivated for making their associated constraints in (|12l) 
tight. In addition, a component in C\ can become deactivated due to a constraint in (|12[) only if it has already lost all 
its children, i.e., each of these vertices in the component is already connected to d 2 . Therefore, by the construction of 
F 2 , each of the marked vertices not spanned by F{ must be connected to di and spanned by F 2 . Hence, the algorithm 
produces a feasible, heterogeneous spanning forest. 

Also, because of the way in which we prune the edges in the final step of the algorithm, if a vertex with label C is 
not spanned by F{, then any other vertex with label C C C is also not spanned by F{. Let X denote the set of vertices 
not spanned by F[. Based on the label of each vertex in X, X can be partitioned into disjoint, deactivated components 
Ci, C2, • • • , C m where each Cj denotes the maximal label of its respective component. 

Consider any deactivated component C I. d can get deactivated during an iteration only if Ci does not have 
children and J2scc- ^i(^) = w(d) — Bound(d) — J2scc- ^(S). Note that d could have lost all its children only 
if all the targets in d are already connected to d-z in F%. Also, during the iteration when C, gets deactivated, no 
target u <E d is connected to any other target v £ T\ Ci in Fi. From corollary HJ it follows that for any u £ Ci, any 
adjacent vertex of u in the deactivated component of C 2 can either be d 2 or another vertex v S Ci. Therefore, during 
the construction of F 2 all the edges that are incident on any vertex already spanned by F[ can be dropped. Hence, any 
vertex spanned by the edges in F[ is not spanned by the edges in F 2 and vice versa. ■ 

The main result of this article is in the following theorem. 

Theorem III.l: The primal-dual algorithm produces a tree with edges denoted by F[ for the first vehicle and a tree 
with edges denoted by F 2 for the second vehicle such that the cost of the edges in these trees is bounded by the cost for 
the dual problem, i.e., 

™ s ti + cost l < 2 E 

e£F[ eEF^ SCT 

Since 2 J^sct ^i(^) ^ s a l° wer bound to the optimal cost of the 2DHTSP, it follows that the cost of the HSF found by 
the primal dual algorithm is at most equal to the optimal cost of the 2DHTSP. 

B. Proof of the Approximation Ratio 

As mentioned in lemma [2] let X denote the set of vertices not spanned by F[. Based on the label of each vertex in X, 
X can be partitioned into disjoint, deactivated components C\, Ci, ■ ■ ■ , C m where each Ci denotes the maximal label of 
its respective component. In order to prove the above theorem, we first simplify the dual cost obtained by the algorithm 
as follows: 



2 y £Y 1 (S) = 2 J2 Y 1 (S) + 2J2J2 Y ^ 

SQT SCT,S<£Ci,i=l,..,m i=1 SQCi 

m 

= 2 £ Y 1 (S) + 2J2 Y,Y2(S). (15) 

SC.T,SgO it i=l,..,m i=1 SCCi 

Now, we express the cost of the edges in the first tree in terms of the dual variables as follows. Note that edge e is added 
to Fi and consequently appears in F{ only if the corresponding constraint in (1101) is tight, i.e., costl = X)s ee5i(S) Yi(S). 
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Therefore, 



E cos< ' = E E Y ^ s ) 

= ^y 1 (s-)Kf|^(5')l- 

SCT 

Since F{ p| <5i(5) = for any S C C;, we can further simplify the above equation to 

Y, co <= E ^M'pi^os)!. (i6) 

e&Fl SCT,S<£Ui,i=l,..,m 

Similarly, we can also express the cost of the edges in the second tree in terms of the dual variables as follows. From 
lemma [2] note that F' 2 can be decomposed into a set of disjoint sets F' 2i where each consists of edges that form a 
tree spanning each target from Ci and the depot d 2 . An edge e is added to F2 and consequently appears in F 2i only if 
the corresponding constraint in pip is tight, cost 2 e = ^2g. e£ s 2 {S) sec - ^(S) where d2i(S) consists of all the edges with 
one endpoint in S and another end point in Ci 1J{g?2} \ S. 



E cost ^ = E E cost e 

m 

= EE E 

i=l eeF^ S:ee5 2l (S),SCCi 
m 

=e e ^(5) i^n^^i. a?) 

i=1 SCCi 

Therefore, from equations (fl5l [T6l fT7| . the proof for the theorem reduces to showing the following result: 

m 

y, yiCfioi^n^^i+E e y ^ i^n^(5)i as) 

SCT,S£Ci,i=l,..,m i=1 SCC; 

m 

<2 £ ^i(5) + 2^ £ y a (S). (19) 

SCT,S^C„i=L..,m i=1 SCCi 

The above result can be shown by proving that during any iteration, the increase in the primal cost (the left-hand side 
of the above inequality) is at most equal to the increase in the dual cost (the right-hand side of the above inequality). 
To see this, let us choose any iteration of the primal-dual algorithm. At the start of this iteration, let N a be the set of all 
the active components in C\ such that each active component in this set is not a subset of X and Nd be the set of all the 
inactive components in C\ such that each inactive component in this set is not a subset of 3C. Note that one of inactive 
components of Nd must consist of the depot d\. For i = 1, • • • , m, let M a i denote the set of all the active components 
in C2 such that each active component in this set is a subset of CV Also, let Md denote the inactive component in C 2 
that consists of the depot di- 

Now, form a graph Hi with components in N a {jNd as its vertices and edges e G F{f]Si(C) for C G N a \jNd as 
edges of Hi. Similarly, form a graph H 2 i with components in M ai [J Md as its vertices and edges e G F^f] 62(C) for 
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C G M a i\J Md as edges of H 2 i. Let deg(y,G) represent the degree of vertex v in graph G. During the iteration, the 
dual variable corresponding to each of the active components is increased by £ m in- As the result, the left hand side of 
the inequality will increase by £ m in(Ylv&N deg{v, Hi) + YlT=i Ylv&M deg(v, i?2i)) whereas the right hand side of the 
inequality will increase by 2e m i n (N a + Y^iLi M a i)- Therefore, basically, the proof is complete if we can show that 

m m 

deg(v, H 1 )+J2 12 de ^ v ' H ^ - 2 ( N * + 12 Mai ^ ( 20 ) 

v£N a i=l vGM ai i=l 

We now claim that any vertex v in H± that represents an inactive component in Nd must have its degree deg(v, Hi) > 2 
unless the inactive component contains the depot d±. This result follows from the fact that a component, which does 
not contain di, can become inactive in Ci only if the constraint associated with this component in (|12[) becomes tight. 
Therefore, all the vertices in this inactive component must be marked. Also, if vertex v is a leaf (deg(v, H\) = 1) then 
pruning all the edges from this inactive component will not disconnect any unmarked target from d±. Hence, the pruning 
step of the algorithm will ensure that an inactive component can never be a leaf vertex in Hi unless it contains d±. Now, 
we show the final part of the proof: 

m 

J2 deg(v,H 1 ) + J2 12 deg(v,H 2i ) (21) 

veN a i=l v£M ai 

= Yl deg(v,Hi)- Y degiv,^) (22) 

v&N a U N d v£N d 
m 

+ £[ 12 deg(v,H 2l )-deg(M d ,H 2t )} (23) 

»=1 v£M ai \J{M d } 

rn, 

<2(\N a \ + \N d \ - 1) - (2\N d \ - 1) + 2 £ \M ai \ (24) 

i=i 

m 

<2|JV |+2^|M oi |. (25) 
i=i 

Hence proved. 

IV. Approximation Algorithm for the 2DHTSP 
The approximation algorithm for the 2DHTSP is shown in algorithm |U 
Algorithm 4 : Approximation algorithm for 2DHTSP 

l: Compute a heterogeneous spanning forest using the primal-dual algorithm. Let the two trees corresponding to the 

first and the second vehicle be denoted as Treei and Tree 2 respectively. 
2: For i — 1,2, do the following: 

• Form an Eulerian graph, for the i th vehicle by doubling the edges in Treei. 

• Using €i, find a walk that visits each of the edges in <Bi exactly once. Next, shortcut this walk to find a tour for the 
jth ve hi c i e suc h that each vertex in is visited exactly once. 



Corollary 2: Algorithm @] presented for 2DHTSP has an approximation ratio of 2. 
Proof: Since the costs satisfy the triangle inequality, the sum of the cost of the tours obtained from the algorithm 
is at most equal to twice the cost of the edges in the trees. Therefore, from theorem IIII.ll it follows that the cost of 
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the feasible solution obtained is at most equal to twice the optimal cost of the 2DHTSP. In addition, the computational 
complexity of the algorithmic is dominated by the primal-dual procedure which can be implemented in |T| 2 log \T\ steps. 
Therefore, the approximation ratio of algorithm |4] is 2. Hence proved. ■ 
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